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ABSTRACT 

We present five band imaging of the Vega debris disc obtained using the Herschel Space Observatory. These data span a wavelength range of 
70-500 yum with full-width half-maximum angular resolutions of 5.6-36.9". The disc is well resolved in all bands, with the ring structure visible 
at 70 and 160 /jm. Radial profiles of the disc surface brightness are produced, and a disc radius of 11" (~ 85 AU) is determined. The disc is seen 
to have a smooth structure thoughout the entire wavelength range, suggesting that the disc is in a steady state, rather than being an ephemeral 
structure caused by the recent collision of two large planetesimals. 
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1. Introduction 

Debris discs, of which the Vega (a Lyrae) disc is the archetype, 
] are characterised as discs of dusty material generated by the col- 

■ Usion of planetesimals in belts surrounding main sequence stars. 
' The ages of the s tars which exh ibit these discs (~ 350 Myr in the 
, case of Vega; Song et al.ll2000h precludes the possibility for this 

dust to b e primordial, as the time scale to remov e such dust is 
, <10Mvr (iBackman & Paresc3ll993HWvatil2008h . 

The debris disc around Vega was first detected by 
' lAumann et al.l (1 19841) as an infrared exce ss using the Infrared 
I Astronomical Satellite (IRAS: iNeugebaue r et al.lll984l) . and has 
- been extensively studied in th e infrared and subm i llimetre over 
' the subsequent 25 years (e.g. [Holland etal.iri998l: IWilner et al.l 
! I2002t ISu et al.ll2005t iMarsh et al.ll2006h . The appearance of the 

■ disc has been found to vary significantly across this wavelength 
I regime, changing from a smooth axisymmetric structure in the 
. infrared (ISu et al.ll2005l hereafter S05), to a structure in the sub- 

■ millimetre, w herein the majority o f the emission lies in two dis- 
] Crete clumps (lHollandetal.lll998h . 

In order to understand the reason for the variation in struc- 
' ture with wavelength it is important to first understand the ori- 
gin of the clumps seen in the submillimetre. The recent colli- 
sion of two massive planetesimals is one option, however, given 
the age of Vega, the statistical likelihood of this occurring with 
two bodies of sufficie nt mass to expla in the submillimetre ob- 
servations is low ([Wvatt & Dent 20021) . A mo re favourable al- 
ternative, firs t proposed bv lWilner et al.l(l2002h and m odelled by 
IWvattI (120061 hereafter W06) and lReche et al.l (l2008h . is that the 
clumps are dust grains trapped in resonance with a planet near to 
the disc. In this scenario the large dust grains (larger than a few 
mm) are trapped in these resonances, while smaller intermediate 



sized grains (a few//m-mm), having been perturbed by radiation 

pressure, have a more uniform distribution in the diso 

Recent analysis and modelling of Spitzer JWerner et al.l 

l2004h mid-infrared data have reached contradictory conclusions. 
S05 find the disc to be ephemeral; in this scenario the disc is the 
result of a recent massive collision of planetesimals, and the sub- 
sequent collisional cascade. This results in a high mass of very 
small grains (less than a few pm) which are blown out of the 
system by radiation pressure immediately upon creation, result- 
ing in the large disc extent observed. Conversely, Miiller et al] 
(I2OIOI) succeed in reproducing the surface brightness radial pro- 
file using intermediate size grains in elliptical orbits around the 
parent planetesimal ring, and therefore conclude that it is con- 
sistent with a steady-state model. In the steady-state model, dust 
that is destroyed, either by being drawn in to the star due to 
Foynting-Robertson drag or blown out of the disc by radiation 
pressure, is continuously replenished by a steady collisional cas- 
cade within the planetesimal belt. 

If the small blown-out grains are the origin of the emission 
observed in the mid-infrared then W06 predicts that spiral fea- 
tures, emanating from the submillimetre clumps, should be visi- 
ble with high-resolution imaging; a smooth structure would sup- 
port the steady-state model. 

In this paper we present five-band far-infrare d imaging of 
the V ega debris disc obtained with the Herschel dPilbratt et al.l 
2010|) Photodetector Array Camera and Spectrometer (PACS 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



Poglitsch et al . 7010) and Spectral and Photometric Imaging 
Receiver (SPIRE; Griffin et al. 201®. We discuss the initial anal- 
ysis and disc parameterisation, and relate these results to the 
ephemeral and steady-state disc models. In Section 2 we present 
the Herschel data, outline the processing performed, and analyse 
the disc structure and properties. These data are then compared 
with results from the recent Spitzer obser vations (S05) and disc 
modelling of W06 and lMiiller et al.l(l2010h . with our conclusions 
summarised in Section 3. 
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Fig. 1. Data for the Vega debris disc from 70//m to 500 jum from left to right respectively (top row - star-subtracted images; bottom 
row - modelled images). All images are scaled linearly, and both images within a given band are equally scaled. The white circle 
represents beam FWHM in each band, and the contour lines are in steps of 5% of the peak flux. The black dashed circle represents 
the location of the disc at a radius of 11". 



2. Observations and data processing 

We obtained images of Vega and its associated debris disc at 
70 and I60fim with PACS, and 250, 350, and SOOfim with 
SPIRE. The data cover an angular scale of ~25 and 64 sq. ar- 
cmin. for PACS and SPIRE respectively, with beam full-width 
half-maxima (FWHM) of -5.6, 11.3, 18.1, 25.2, 36.9" for the 
short to long wavelength bands. The data were obtained in scan- 
map mode for both instruments using the nominal observing pa- 
rameters in both cases, and scanning rates of 10 and 30" per 
second for PACS and SPIRE respectively. The total on-sky ob- 
serving time was 5506 s and 6120s and comprised 70 and 16 
map repetitions for PACS and SPIRE respectively. 

The PACS data were high-pass filtered to remove low- 
frequency noise using a cut-off scale of 3.7'. The data were then 
binned to a map using the default photProject met hod in th e 
Herschel interactive processing environment (HIPE: IOtBl2010l) . 
The pixelisation of these maps was set to 1 and 2" per pixel, 
equivalent to ~l/5 of a beam. The maps have 1-cr noise levels of 
0.045 and 0.083 mJy arcsec"-, which includes 10 and 20% flux 
calibration uncertainties, in the 70 and 160 yum bands. The back- 
ground was removed from the maps by subtracting the median 
in the local vicinity of the source. 

The SPIRE data were also reduced using HIPE and maps 
were obtained via the default na'iveMapper task. The 16 re- 
peat observations allowed the data to be binned to 4, 6, and 
9" sized pixels without losing complete sampling across the 
source. The 1-cr noise level obtained was 0.014, 1.8 x 10 "*, and 
2.8 X 10 "^mJy arcsec - in the 250, 350, and 500 Afm bands re - 
spectively; the calibration error is ~15% JSwinvard et alj|2010l) . 
The SPIRE beams exhibit a 1.07, 1.12, and 1.09 ellipticity, as 
described in the SPIRE beam model release note. 

In order to more easily assess the structure in the disc, the 
photospheric contribution from the star was subtracted from the 
image (top row Fig. [T]i- This was achieved by scaling a high 
signal-to-noise ratio observation of Vesta, which can be regarded 
as a point-source image, to an appropriate flux level to obtain a 



model for the stellar contribution in each band. The Vesta image 
was rotated to match the position angle of the telescope used in 
the observation of Vega before subtraction of the model star The 
ph otospheric flux in each band was estimated using data given 
bv lRieke et aP (l2008h . and colour corrections of 1 .02, 1 .07, 0.96, 
0.99 a nd 1.04 were applied f or the 70 to 500 Afm bands respec - 
tively jPoditschetal'lboid: ISPIRE Observers' Manual |20 Id) . 
The resulting flux density estimates were 793, 162, 63, 32 and 
16 mJy for the 70 to 500 jum bands respectively. The stellar mod- 
els were then subtracted from the reduced data. 

The beam model used in star subtraction for both the PACS 
and SPIRE data was obtained using the same observing mode 
parameters as the original data. The PACS beam exhibits a char- 
acteristic tri-lobe structure, while the SPIRE beam contains clear 
side-lobe structure and a small ellipticity. Uncertainties in these 
beam models result in star subtraction artifacts in the central re- 
gion of the image. This is especially pronounced in the 70jum 
image where the stellar contribution is highest. Consequently the 
star-subtracted maps can only be used to assess the disc structure 
at radii larger than 5 arcsec (~40 AU). 

2.1. image analysis 

The five star-subtracted maps are presented in the top row of 
Fig. [U and range from 70-500 /^m from left to right. The disc 
is resolved in all bands, and shows a smooth and axisym- 
metric structure; centroids of the stellar and disc components 
show less than 1 pixel difference in the position of the star 
with respect to the disc in the 70 and 160yum bands, where 
the star and disc locations can be readily identified. The 30% 
contour of the five star-subtracted discs were fit by an ellipse 
and were found to be extremely circular, with ellipticities of 
~1.01±0.002, 1.02+0.003, 1.04±0.04, 1.03+0.06 and 1.11+0.09 
for 70-500 yum. The 500 //m image exhibits a significantly higher 
ellipticity than the other bands, however, the relatively low res- 
olution in this band coupled with the 1.09 intrinsic beam el- 
lipticity make the statistical significance of this measurement 
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too low to draw a robust conclusion. The flux densities, ob- 
tained via aperture photometry, for the star plus disc system, are 
10.12+1.18, 4.61+0.9, 1.68+0.26, 0.61+0.10 and 0.21+0.04 Jy 
from 70 to 500 //m respectively, with the calibration error domi- 
nating the uncertainty. These measurements agree well with in- 
tegrated ineasurementsjnadeat_sim^^^ by other fa- 
cihties (iMarsh et al.ll2006t ISu et al.ll2006l) . A single aperture of 
radius 30" was used for all bands, with the same aperture ran- 
domly placed around the source to quantify the noise. 

With the star-subtracted from the image the ring structure of 
the disc can be clearly identified at 70 and 160/im. The disc is at 
a radius of ~H", which at the distance of Vega (7.76 pc) corre- 
sponds to ~85 AU. This equates to ~4 beam half width half max- 
ima from the central star, making this detection robust against 
artifacts from the star subtraction. This is in agree ment with pre- 
vious infrared (S05), and submillimetre estimates ( iHolland et alj 
Il998t iMarsh et aDl2006l) . The inner cavity is not visible in the 
SPIRE data due to the decreased resolution relative to the shorter 
wavelength bands; the large scale disc size remains comparable. 

2.2. Radial profiles and surface brightness modelling 

The face-on nature of the Vega disc allows us to obtain data on 
the general disc structure and extent by azimuthally averaging 
the radial intensity profiles. 805 performed such an analysis for 
the 24, 70, and 160//m data from Spitzer and find that the disc 
profile can be fitted by r^^ and r "^ power laws for the inner and 
outer disc respectively. Radial profiles for Vega, derived from 
the higher resolution Herschel data, are presented in Fig. 2 for 
the raw and star-subtracted maps. The stellar model used for star 
subtraction is plotted for reference. Radial step sizes equal to the 
map pixel scale were used out to a radius of 90". 

The disc extends to a radius of ~r in all bands before the 
signal-to-noise ratio drops and the data become subject to un- 
certainties in the baseline removal. The disc radius is in agree- 
ment with that found by S05 in both 70 and 160/im bands, 
with differences 6.4 and 2.8% respectively. To compare these 
data the PACS data were convolved to the Spitzer resolution 
and the disc radius was measured at a surface brightness of 
0.5 mJy arcsec". The PACS radii measured at the Spitzer resolu- 
tion were 34.7++j}^ " and 31.6!i^| " for the 70 and 160 /^m bands 
respectively, with the errors based on the flux calibration uncer- 
tainties. 

The drop-off^ in radial profile at high radius appears linear in 
the log-linear plots in Fig. 2. The functional form of this slope, 
characterised between radii of 20-50", is log[Q(5'v) = -0.63r + 
1.61 at lOjjm, where Sv is, the surface-brightness at radius r. 
This is in contrast to the power-law slope identified by S05 for 
the PACS bands. However, S05 fit the disc out to larger radii, 
and within the region we fit the same functional form could be 
similarly applicable. An accurate comparison with S05 at radii 
larger than ~50" if difficult as these data are highly affected by 
uncertainties in the background subtraction. 

The clearly defined ring identified in Fig. [1] is evident again 
in the radial profiles at 70 and 160jum, and is defined by the 
peak and turn-over of the disc profile at a radius of 11". The 
structure observed, however, makes it difliicult to obtain a dis- 
crete measurement of the inner and outer edges of the disc. As 
an alternative characterisation we measured a half-width half- 
maximum (HWHM) size of the disc, outward in radius from the 
peak disc brightness at 11" (Fig. 2), and using this peak as the 
reference maximum. We obtained HWHM sizes for the 70 and 
160 fim bands, in which the disc radius is identifiable, of 9 and 
11.3" respectively. 
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Fig. 2. Radial Surface brightness profiles in the five Herschel 
bands for the original, the star-subtracted, and modelled data, 
as well as the scaled stellar photospheric model. Error bars of 
1-cr are given for the star plus disc measurements. 

To further assess the potential disc structure, we developed 
a simple model based only on the observed surface brightness 
properties at 70jum. The model comprises two different surface 
brightness distributions, separated at a reference radius, tq. The 
inner profile, r < ro was chosen to be a Gaussian function peak- 
ing at ro, and the outer profile, ro < r, was defined using the 
functional form found for the 70jum outer disc described above. 
The distinction between the inner and outer discs, as parame- 
terised in this model by ro, is a purely observational definition, 
based on the transition from one brightness distribution model to 
the other. A more physically meaningful definition for the sep- 
aration between the inner and outer disc is the peak of the disc 
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Fig. 3. Difference image for the star-subtracted data and the 
model image at IQfim. In this instance the model image is noise- 
less. 

profile found at a radius of 11" as described above. The model 
was then convolved with the appropriate beam model and com- 
pared to the measured surface brightness profile. 

A best-fit model was found with parameters ro = 14" and 
inner Gaussian FWHM = 20". Models for the bands longward 
of 70yum were created using the same intrinsic surface bright- 
ness profile and convolved with the appropriate beam. The peak 
value of the resultant models were scaled by 0.33, 0. 10, 0.04 and 
0.01 for the 160 to 500 //m bands respectively to match the data, 
and placed in an empty region of the original map to replicate 
the realistic instrumental noise (bottom row Fig. [Til. The radial 
profile of the output modelled image was measured and plotted 
with the real data in Fig. 2. There is good agreement between 
the modelled and real radial profiles, with l-cr residuals in the 
20-50" fitted region below 0.05 mJy per sq. arcsec across all five 
bands, implying that the underlying structure of the disc across 
the wavelength range is similar. 

As the model data are known to be perfectly smooth, a direct 
comparison can be made with the real data to assess the signifi- 
cance of any potential small scale structure seen within the disc. 
For example, at 160yum very low-level structure can be seen, 
with flux enhancements in the northern and southern parts of the 
ring. However, the counterpart model disc image in Fig.[T]shows 
similar features, indicating that these are at the level of the noise, 
and should not be attributed to true disc structure. 

Subtraction of this simple uniform model, without added 
noise, from the original data provides another method to easily 
identify disc structure which is otherwise difficult to detect in the 
presence of the larger disc; a difference map at lO/im, the high- 
est resolution band, is shown in Fig. |3] The brightest structure 
in this image lies within the inner disc, where the beam subtrac- 
tion artifacts are strongest. This can also be seen in Fig. 2 as a 
difference between the radial profiles of the model and the star- 
subtracted image at ~6". The dark features seen in the outer disc 
region also correspond in structure and position to the triple- 
lobed beam pattern, and therefore are disregarded as potential 
disc features. With the exception of the beam subtraction struc- 
tures, there is no sign of any other clumpy structure associated 
with the disc, down to the noise limit of these data. 

3. Discussion and conclusions 

The structure observed in the Herschel data shows no sign of 
clumps. There are also no visible spiral arm features, predicted 
by W06, if the disc emission at mid-to far-infrared wavelengths 
is dominated by small blown-out dust grains (W06 Fig. 3; right 
hand panels, /? =1-10). The smooth structure observed is most 
consistent with the steady-state model, wherein the emission is 



dominated by intermediate size dust grains in elliptical orbits 
about the parent planetesimal belt. This model was found to 
simultaneously give good agreement to the data in all bands, 
which is unexpected, as the more distant grains should have a 
lower temperature, and suggests that the mean grain size de- 
creases with distance from the star. This is in-keeping with 
observational data which shows a larger disc at shorter wave- 
lengths. Full modelling of the radial grain size distribution will 
be presented in Sibthorpe (2010, in prep). 

1 . We presented resolved images of the Vega debris disc system 
in five bands ranging from 70-500 jt/m obtained using the 
Herschel PACS and SPIRE instruments. 

2. The peak surface brightness of the dust disc was identified at 
70 and 160 A^m at a radius of 1 1" (85 AU). 

3. The surface brightness profile was found to be well fit in the 
outer disc by a log[Q(5'v) cc -0.63r distribution, with a differ- 
ent scale factor at each band. The inner profile (r < ro) was 
likewise modelled, with a Gaussian profile of FWHM = 20" 
found to provide a good fit. The change in surface bright- 
ness distribution, occuring at a radius of ~14" (~109 AU), is 
used to observationally define the distinction between the in- 
ner and outer disc. This model was found to simultaneously 
give good agreement to the data in all bands. 

4. The structure of the disc was found to be smooth, with no 
clumpy structure to the sensitivity limit of these data. 

5. While these data cannot preclude the option that the Vega 
disc is the result of a large planetesimal collision, making it 
ephemeral in nature, these data support the hypothesis that 
the Vega disc is steady-state in nature. 
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